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Recent X-ray observations by Fermi/GBM discovered a new torque reversal of 4U 1626— 67 after 18 years of 
steady spinning down. Using Swift/BAT observations we were able to center this new torque reversal on Feb 4 
2008, lasting approximately 150 days. From 2004 up to the end of 2007, the spindown rate averaged at a mean 
rate of ~ v = —4.8 X 10 -13 Hzs -1 until the torque reversal reported here. Since then it has been following 
a steady spinup at a mean rate of ~ u = 4 X 10 -13 Hzs -1 . The properties of this torque reversal, as well as 
the lack of correlation between the X-ray flux and the torque applied to the neutron star before this transition, 
challenges our understanding of the physical mechanisms operating in this system. 



1. INTRODUCTION 

The accreting-powered pulsar 4U 1626-67 was dis- 
covered by Uhuru [l7j . This low mass X-ray binary 
(LMXB) consists of a 7.66 s X-ray pulsar accreting 
from an extremely low mass companion (0.04 M© 
for i — 18°) [30( | . Although orbital motion has never 
been detected in the X-ray data, pulsed optical emis- 
sion re pro cessed on the surface of the secondary re- 
vealed [36( the 42 min orbital period. The faint op- 
tical counterpart (KZ TrA, V~17.5) has a strong UV 
excess and high optical pulse fraction [H, HH • A per- 
sistent 48 mHz quasi-periodic oscillation (QPO) has 
been detected in the X-ray emission [24|, |46| • [I?! in- 
ferred a neutron star magnetic field in the range (2.4- 
6.3) xlO 12 G. To compute this magnetic field range 
a source distance of 5-13 kpc was assumed @, and 
references therein]. A ^37 keV absorption cyclotron 
feature was found in the 0.1-200 keV BeppoSAX spec- 
trum [37J . 

For more than a decade after the discovery of pul- 
sations [44| the source underwent steady spin-up at a 
mean rate of ~ v = 8.5 x lO^Hzs" 1 [1] (see FigQ] 
Top). Monitoring of the source by the Burst and 
Transient Source Experiment (BATSE) on board the 
Compton Gamma Ray Observatory (CGRO) starting 
in April 1991, found the pulsar spinning down, im- 
plying a changed sign in the accretion torque 0, 51] ■ 
During the 7 years after the first torque reversal, the 
pulsar spun-down at a rate of ~ v = —7.2 x 10~ 13 Hz 
s" 1 i. 

We present a long term timing and spectral analy- 
sis using all the available Fermi/GBM data since its 
launch in 2008 June 11 and over 5 yr of hard X-ray 
Swift/BAT data from2004 up to 2009. 



2. Fermi/GBM 

2.1. OBSERVATIONS 

Since 2008 June 11 4U 1626-67 has been contin- 
uously monitored by the Gamma-ray Burst Monitor 
(GBM)(35|, on board the Fermi observatory. Tim- 
ing analysis was carried out with GBM CTIME data, 
with 8 channel spectra every 0.256 seconds. The total 
exposure time was ^13.75 Ms. The GBM is an all- 
sky instrument sensitive to X-rays and gamma rays 
with energies between ~8 keV and ~40 MeV. GBM 
includes 12 Sodium Iodide (Nal) scintillation detec- 
tors and 2 Bismuth Germanate (BGO) scintillation 
detectors. The Nal detectors cover the lower part of 
the energy range, from 8 keV to about 1 MeV. The 
BGO detectors cover the energy range of ~150 keV 
to ^40 MeV. Only data from the Nal detectors were 
used in the analysis. 



2.2. TIMING ANALYSIS AND RESULTS 



All intervals of CTIME data from the 12 Nal 
detectors are selected for analysis after excluding 
those containing high voltage transients, phosphores- 
cence events, rapid spacecraft slews, South Atlantic 
Anomaly induced transients, electron precipitation 
events and gamma-ray bursts (see Mark H. Finger et 
al., these proceedings). Source pulses are then sepa- 
rated from the background by fitting the rates in all 
detectors with a background model, and subtracting 
the best fit model. Then we combine the residuals 
over detectors with time dependent weights which are 
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Figure 1: (Top). Pulse frequency history of 4U 1626-67. 
(Bottom Left). Zoom of the top panel, with Fermi GBM 
pulse frequency measurements since 2008 August. A 
change in the sign of the torque was found after 18 years 
of the source spinning down. (Bottom Right). Zoom of 
the top panel, with SwiftfBAT pulse frequency history 
covering the 2008 torque reversal. Error bars are smaller 
than the plotted symbols. 
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Figure 2: Top panel. Swift/BAT spin-up rate history of 
4U 1626-67 . Middle panel. Average 15-50 keV BAT 
count rate vs. time. Error bars are smaller than the 
plotted symbols. Bottom panel. BAT count rate vs. 
spin-up rate for all the period (circles). A correlation 
pattern is observed specially during the torque reversal 
(only square symbols). 



3. Swift/BAT 

3.1. OBSERVATIONS 



proportional to the predicted (phase averaged) count 
rates from the pulsar. 

Short intervals (~300s) of these combined residuals 
are then fit with a constant plus a Fourier expansion to 
determine a pulse profile. The profiles are divided into 
six day intervals and the pulse frequency and mean 
profile determined in each interval with a search of 
pulse frequency for the maximum of the Y n (n=2) 
statistic |13j . 

Our monitoring of 4U 1626-67 with Fermi/GBM 
starting in 2008 August, discovered the pulsar 
spinning-up [gj . Fig Q] (Right) shows the pulse 
frequency history using data from this monitoring. 
4U 1626-67 seems to be increasing in v. Follow-up 
Fermi/ GBM observations confirm that the pulsar it is 
currently spinning-up at a mean rate of ~ v—Ax 10~ 13 
Hzs" 1 . 



The Swift Gamma-ray mission [15l | was launched 
on 2004 November 20. The hard X-ray (15-150 keV) 
Burst Alert Telescope (BAT) on board Swift, monitors 
the entire sky and it produces continuous streams of 
rate data. For the present study we have analyzed 
more than 4 years of BAT quadrant rates observa- 
tions (1.6 sec sampling; four energy bands; four sepa- 
rate spatial quadrants) when 4U 1626-67 was visible 
(total exposure time ~ 13 Ms). For the hardness ra- 
tio analysis we used count rates from the Swift/HAT 
transient monitor results provided by the Swift/BAT 
team 1 . 

3.2. TIMING ANALYSIS AND RESULTS 

A similar procedure was followed for the Swift/BAT 
quadrant rates timing analysis. Initial good time inter- 
val (GTI) files are obtained using the maketime ftool 



1 http: / /hcasarc. gsfc.nasa.gov/docs / swift / results /transients 
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(heasoft-6.6.1) 2 . Then a filtered version of the quad- 
rant rates is obtained, to then finally be barycentered 
using the ftool barycorr. With the ftool batmasktgimg 
the pixel exposure fraction for each quadrant is com- 
puted for the center of each (refined) GTI interval. 
Pulse profiles for each good GTI interval are com- 
puted. First the rates for each quadrant are fit to a 
quadratic+Fourier expansion. Then the Fourier co- 
efficients are combined using the quadrant exposures 
to produce mean profiles (with units of counts s _1 
cm -2 ). In a final stage, the Y„ (n=2) statistic is again 
used in intervals of 35 days and a frequency search for 
pulsations is carried out. The spin rates were com- 
puted by fitting a linear function to the frequencies, 
which were divided into 21 time intervals. 

Swift/BAT observations allowed us to cover the 
evolution of this second torque reversal. We found 
that the pulsar spun-down at a mean rate of ~ i> = 
—4.8 x 10 _13 Hzs -1 until the source reversed torque. 
Fig [1] (Left) shows that the transition took place at 
around MJD 54500 (2008 Feb 04) and lasted approx- 
imately 150 days. In the bottom panel of Fig [3] we 
can see that there is a strong correlation between the 
Swift/BAT count rate and the spin-up rate especially 
during the reversal. We have not observed any sig- 
nificant change in pulse shape, not even during the 
reversal. 



4. RXTE 

4.1. OBSERVATIONS 

The Rossi X-ray Timing Explorer {RXTE) [| 
carries 3 instruments on board. The Proportional 
Counter Array (PC A) [p (2-60 keV), the High En- 
ergy X-ray Timing Experiment (HEXTE) fl8| (u p to 
200 keV) and the All Sky Monitor (ASM) [2g (2- 
10 keV). Two RXTE/PCA observations from 2008 
March 5 and 13 were used (ID 93431-01-01-00 and 
93431-01-02-00; 7.174 ksec). For spectral analysis 
we selected PCA Standard-2 data and HEXTE Stan- 
dard Modes (Archive) data (64-bin spectra every 16s). 
For the long-term hardness ratio analysis we used the 
ASM daily flux averages in the 1.5-12 keV energy 
range from the HE AS ARC archive 3 . 



4.2. SPECTRAL ANALYSIS AND 
RESULTS 

RXTE/PCA (2.5-20 keV) and HEXTE (18-100 
keV) spectra were fitted in XSPEC 11.3.2 with 
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2 http: / /heasarc. gsfc.nasa.gov/docs/software. html 
3 http: / /heasarc. gsfc.nasa.gov/docs/archive. html 



Figure 3: Hardness-intensity diagram of 4U 1626-67. For 
HR, Swift/~BAT count rates (15-50 keV) were selected as 
the hard band and RXTE/ ASM count rate (1.5-12 keV) 
as the soft band. During the reversal a transition from 
hard to soft is seen (green circles). 



two models used by [41| . Using these models al- 
lows us to compare our spectral study with previ- 
ous works by [H, H^, [3?], El and update the long- 
term X-ray flux history of 4U 1626-67 relative to 
the flux measured by HE AO 1 (Chakrabarty et al. 
(1997); Krauss et al. (2007)) . The first model in- 
cludes a low-energy absorption, a blackbody compo- 
nent, a power law and a high-energy cutoff at ~20 
keV (WABS (GAUSS+BBODY+POWLAW) HIGH- 
ECUT). A broad line near 6.5 keV significantly im- 
proves the present fit and indicates the presence of an 
iron line, also suggested by [4l| in their ( 0.7-100 keV) 
spectral analysis of this source. The column density 
of cool material in the line of sight was fixed in our 
study since it could not be constrained. A value of 
1.3 x 10 21 cm~ 2 was selected from (25[. The spectral 
parameters obtained are shown in Table 1. We fit 
in addition the same model with a bremsstrahlung in- 
stead of a blackbody component , obtaining a compat- 
ible fit. Table U summarizes the spectral parameters 
obtained. 



5. HARDNESS RATIO ANALYSIS 

Fig [3] shows Hardness-intensity diagram (HID) for 
4U 1626-67 using BAT count rate. The hardness ra- 
tio (HR) was defined as the ratio 15-50keV/1.5-12keV 
(BAT/ ASM). To reduce large uncertainties the light 
curves were rebinned and then the HR were computed. 
This allow us to study the long-term spectral variabil- 
ity of 4U 1626-67 , including the transition, since the 
2 RXTE observations do not provide us any direct 
comparison between before and after the torque re- 
versal. From that figure we can see that there is a 
transition from hard to soft during this new reversal 
of 4U 1626-67. 
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Table I RXTE/PCA and HEXTE CONTINUUM SPECTRAL FITS" 



Observation 
(MJD) 



^cut 



El 



Gaussian 



Gauss. 



Gauss. 



T t 

- 1 - BBody 
T t 

Bremss 



norm Flux"* x. 2 (DOF) 



B Bod 



93431-01-01-00 0.75(2) 18.2 
(54530) 0.74(2) 18.2 



+0.1 
-0.3 
+0.3 
-0.2 



8.5(4) 6.05 
8.5+0 4 6.42 



+0.5 
-0.17 
+0.17 

-0.4 



1.6+°; 3 2 2.4 



1.4 



+0.1 

-0.2 



1.5 



+ 4 
-0.8 
hi. 2 
-0.6 



0.615+o°°| 0.0013 1.01(8) 1.15(114) 
1.74+°^ 0.108 1.01(1) 1.15(114) 



,+0.19 



93431-01-02-00 0.71 Z Q0 4 17 - 90 ~- o ; 



-1+0.5 e O+0.6 
■ 4 _ .6 b - Z -0A 



(54538) 



0.71+°;^ 17.96 _„ , 



+°' 14 8.4(6) 6.9(2) 



1.5+0 2 2.23 + 1 g 0.654 +o°o°o 4 9 0.0013 1.006(12) 1.29(114) 
0.5+Q 3 °- 396 -o^ils 2 - 47 -o°o 3 8 °- 068 1 -004(9) 1-28(114) 



mod WABS (GAUSSIAN+ BLACKBODY / BREMSSTRAHLUNG +POW) HIGHECUT (Ah = 1.3 x 10 21 cm" 2 fixed, 
uncertainties at 3<r level) 



* Photon Index 
t keV 

" XlQ- 3 



1 (Lumin/10 39 ergs- 1 )(d/10fcpc)- 2 

2 3.05 x 10~ 16 (47rd 2 ) _1 xEmission measure 
*** xl0~ 9 ergcm~ 2 s~ 1 (2-100 keV) 



6. DISCUSSION 

All previous studies of 4U 1626-67 were focused on 
modeling the spin-up torque applied to the neutron 
star from the accreted material. It was widely believed 
that the spin behavior of the pulsar depended mainly 
on variations of the mass accretion rate onto the stel- 
lar surface and therefore, the rate of mass transfer 
between the system components. A dramatic change 
in the power spectra between the last observation of 
4U 1626-67 during the spin-down phase (2003) [H, 
and observations made soon after the new torque re- 
versal has been recently reported by [13], with the 
35-48 mHz QPO no longer being present, and wide 
shoulders on the pulse fundamental appearing. They 
claimed that the observed behavior of the source can- 
not be a simple case of increased mass transfer rate, 
but is also a change in the accretion flow parameters. 

Analyzing the evolution of the source energy spec- 
trum and possible correlation between the torque and 
X-ray luminosity of the pulsar, [5(| proposed a sce- 
nario in which the torque reversal in 1990 is associ- 
ated with a state transition of the accretion disk to a 
geometrically thick, hot and, possibly, sub-Keplerian 
phase. Following this idea one could associate the 
2008 torque reversal with an inverse transition of the 
disk into its previous geometrically thin Keplerian 
phase. However, the reason for such a transition 
is rather unclear since the level of X-ray flux mea- 
sured before and even after the 2008 reversal is smaller 
than that measured during the reversal in 1990. Fur- 
thermore, both reversals have occurred at almost the 
same timescale (about 150 days), which significantly 
exceeds the dynamical timescale in the hot disk in 
which its transition to the ground state is expected. 
The difficulty to fit the transition timescales observed 
in 4U 1626-67 has also been mentioned by (48[, who 
discussed a possibility to explain the torque reversals 
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Figure 4: The X-ray flux history of 4U 1626-67 relative 
to the flux measured by HE AO 1, in the same energy 
band, from previous works (Chakrabarty et al. (1997): 
circles; Orlandini et al. (1998): triangle; Krauss et al. 
(2007): stars) and two recent RXTE/PCA observations 
(unfilled squares) in the 2-20 keV band. The cross point 
is inferred from PCA flux and the fractional change in 
the Swift/BAT rate. 

in terms of the warped disk transition into a retro- 
grade regime. 

A correlation between the torque applied to the neu- 
tron star in 4U 1626-67 and X-ray flux of the sys- 
tem in the above mentioned models has been adopted 
as one of the basic assumptions. To test the valid- 
ity of this assumption using data derived before 1993 
was rather complicated. This is illustrated in Fig- 
ure[H which shows the 4U 1626-67 X-ray flux history. 
We can see all previous flux measurements and two 
RXTE/PCA recent values from the present work (in 
the 2-20 keV band). These values are relative to the 
flux measured in 1978 by RE AO 1 in the same en- 
ergy band 0, H^, [53, and references therein]. The 
cross point before the 2008 reversal has been inferred 
by scaling the PCA fluxes according to the observed 
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change (2.5 factor) in the Swift/BAT rate, since no 
spectral changes across the transition have been ob- 
served according to the present work. 

As seen from Figure HI the X-ray flux during the 
spin-down phase has decreased by a factor of 2. This 
indicates that the mass accretion rate onto the surface 
of the neutron star, M, and, correspondingly, the spin- 
up torque applied to the star [42| . 

K su = M(GM ns r m )^ 2 , (1) 

during this phase have also decreased by at least the 
same value. If the spin-down torque applied to the 
neutron star during this time were constant one would 
expect the pulsar to brake harder at its fainter state 
close to the end of the spin-down phase. However, ob- 
servations show the situation to be just the opposite. 
The spin-down rate of the neutron star during this 
phase has decreased from \v\ ~ 7x 10~ 13 Hzs _1 [H to 
5 x 10~ 13 Hz s _1 (see Figured]), implying that the pul- 
sar was braking harder at its brighter stage just after 
the torque reversal in 1990. According to the equation 
governing spin evolution of an accreting neutron star, 

l^lv = K su - K sd , (2) 

this means that the spin-down torque, K s( \, during 
the spin-down phase has been decreasing simultane- 
ously with the spin-up torque but at a higher rate 
and, therefore, the pulsar spin evolution during this 
time has been governed mainly by variations of K s & 
rather than K su (here / is the moment of inertia of the 
neutron star). This conclusion seriously challenges the 
possibility of modeling the spin history of 4U 1626-67 
solely in terms of variations of M, and suggests that 
the dramatic increase of X-ray flux observed in 2008 
torque reversal may be a consequence rather than a 
reason for this event. 

With the lack of correlation between the X-ray flux 
and the torque applied to the neutron star, modeling 
of the spin-down torque appears to be the main target 
for theoretical studies of the system. Unfortunately, 
this part of modeling of the magneto-rotational evolu- 
tion of neutron stars remains so far a work in progress. 

Perna et al. (2006) proposed a model where si- 
multaneous accretion from a disk onto the neutron 
star, some material from near the disk-magnetosphere 
boundary is ejected, and either escapes from the sys- 
tem or is recycled back into the accretion disk. Their 
model predicts, however, that the luminosity after a 
spin-down to spin-up torque reversal would be higher 
than the luminosity after a spin-up to spin-down 
torque reversal, which is the opposite of what occurred 
for 4U 1626-67 for this new reversal. Moreover, for 
4U 1626-67 they predicted the full spin-down/spin-up 
cycle would take thousands of years, again inconsis- 
tent with the present observations. 

Finally, the spectral evolution of 4U 1626-67 dur- 
ing the torque reversal differs from that expected in 



models which suggest significant changes of the accre- 
tion flow structure in spin-up/spin-down transitions 
[e.g.][43,l50|. As seen from Figure [31 the spectrum be- 
comes the hardest during the reversal and the value 
of the hardness ratio before and after these events 
does not differ significantly. This indicates that the 
recent torque reversal can be associated with changes 
of physical conditions in the inner part of the disk 
or/and in the region of its interaction with the mag- 
netosphere rather than a significant change of the ac- 
cretion flow geometry. The errors of the observations 
are, however, too large for a justification of particular 
transition model. 



7. CONCLUSIONS 

We report on a discovery of a new spin-down to 
spin-up torque reversal in 4U 1626-67. It occurred 
after about 18 years of the pulsar's steadily spinning 
down and was centered on 2008 Feb 4. The transi- 
tions lasted ~150 days and was accompanied by an 
increase in the Swift/BAT count rate of a 2.5 factor 
(^150%). The pulsar spectrum was harder during the 
torque transition than before or after. A strong corre- 
lation between torque and luminosity is inferred only 
during the transition. The spin-up and spin-down 
rates before and after the transition were almost iden- 
tical (~| v |= 5 x 10~ 13 Hzs _1 ). However, the pulsar 
was braking harder at the beginning of the spin-down 
epoch in 1990 than at its end in 2008. Furthermore, 
the spin-down rate during this epoch was decreasing 
simultaneously with the decreasing of the source X-ray 
luminosity. Finally, the spin-down to spin-up torque 
reversal in 2008 occurred at lower luminosity as the 
spin-up to spin-down torque in 1990. These proper- 
ties cannot be explained with existing models and ap- 
pear to be a clue for further progress in understanding 
the mechanism governing the torque reversals in the 
accretion-powered pulsars. 
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